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ABSTRACT 

Context. Mass loss is an important property in evolution models of massive stars. As up to 90% of the massive stars have a visual 
or spectroscopic companion and many of them exhibit mass exchange, mass-loss rates can be acquired through the period study of 
massive binaries. 

Aims. Using our own photometric observations as well as archival data, we look for variations in orbital periods of seven massive 
eclipsing binary systems in the Cygnus OB2 association and estimate their mass-loss rates and stellar parameters. 

Methods. We use a Bayesian parameter estimation method to simultaneously fit the period and period change to all available data and 
a stellar modelling tool to model the binary parameters from photometric and radial-velocity data. 

Results. Four out of the seven selected binaries show non-zero period change values at two-sigma confidence level. We also report 
for the first time the eclipsing nature of a star MT059. 
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1. Introduction 


The change of a period (P) in a binary system has been com¬ 
monly evaluated by measuring the times of minimum of pri¬ 
mary eclipses and then comparing those times to the calculated 
ephemerides for the eclipses. The period change value is then ac¬ 
quired by fitting a polynomial model through the observed minus 
calculated (O - C) data points. The value of this period change 
can indicate mass loss through stellar winds or mass tra nsfer be¬ 
tween the binary components (ISingh & Chaubevll986h . 

Mass loss in binary systems comes mostly from the mass 
transfer via Roche lobe overflow when the primary star’s photo¬ 
sphere crosses the inner Lagrange point and starts to exchange 
mass with its companion. While the mass is transferred, the cen¬ 
tre of mass of the interacting binary system shifts towards the 
mass-receiving star. Due to the conservation of angular momen¬ 
tum, this mass transfer causes a change in the system’s orbital 
period: the period increases while the mass is transferred to the 
more massive star and vice versa. The period also increases when 
ma ss is lost from the system through stellar winds. 

Is ana et ahl (l2012h studied a sample of massive O stars and 
found that over 70% of them exchange mass with a compan- 
ion. The speckle int erferometric survey of Galactic massive stars 
(iMason et al.l l2009l) expanded upon a decade-old survey and 
noted that 75% of the massive stars in clusters and 59% in 
the field have either a spectroscopic or a visual companion. As 
low inclination or long periods make the binary systems hard 
to detect, these values are only the lower limit and, when cor¬ 
rected to the observational bias, the upper limit for binarity frac- 
tion in bound systems can re ach as high as ~ 90 + 10% 

dKiminki & Kobulnickvl2012l) . 

Massive stars have a relatively short lifetime that is influ¬ 
enced by their mass-loss rates. This mass loss determines the 


evolution of their temperature and luminosity, and also the su¬ 
pernova explosion in the final evolution phases. Massive stars 
provide feedback in the form of ultraviolet (UV) radiation, stellar 
winds, and supernovae which impact the evolution of disc galax¬ 
ies by mixing and driving turbulence of interstellar gas and trig¬ 
gering star formation. They also enrich the interstellar medium 
with heavy elements and drive the chemical evolution of galaxies 
and the metallicity evolution of the universe. Therefore, knowing 
the evolution of massive stars gives us insight into the evolution 
of galaxies. The uncertainties in determining the mass-los s rates 
for m assive stars are still large even in the local Universe (ISmithI 
12014 . 

Constraining orbital parameters for massive binary systems 
is important to test the predictions of theoretical models. In or¬ 
der to determine all the parameters of a single binary system, 
one needs to obtain photometric and spectroscopic data at vari¬ 
ous orbital phases, requiring great observational effort. As it is 
time-consuming to obtain those observations and as the eclips¬ 
ing nature is dependent on the inclination of the system, there 
are not many massive binaries with fully constrained parameters 
known. In this work we present the stellar parameters as well as 
the period and the period change values of seven selected sys¬ 
tems in Cygnus OB2. 

Stellar association Cygnus OB2 is one of the most mas¬ 
sive and richest assoc iations known in our Ga laxy, containing at 
least 60 O-type stars (iNegueruela et al.ll2008h . Determined from 
the stellar ev olution models, the age of the cluster is found to 
be 1-7 Myr dWright et al.ll20L§) . It is dynamically young, sub¬ 
structured and unbound association that is located in the molec¬ 
ular cloud complex C ygnus X and its to tal estimated mass is 
(1.6+ 0.3) X 10^ Mr. dWright et al.ll20l4 . 
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Table 1. References to all the data that is used in our period change 
analysis. 


Star 

RA 

DEC 

References 


Hour 

Deg 



A36 

20:34:59 

-h41;36:17 

1,2,3 

B17 

20:30:27 

-h41:13:25 

1,2 

MT059 

20:31:11 

-h41;31:54 

1,2,4,5 

MT696 

20:34:00 

-h41:17:35 

1,4,6,7 

MT720 

20:34:06 

-h41:08:09 

1,4,5 

Schulte 3 

20:31:38 

-h41:13:21 

1,2,4,5,8 

Schulte 5 

20:32:22 

-h41:18:19 

1,2,4,8,9,10, 




11,12,13,14,15,16 


References. ( 1) iT elescope (our data set 1); (2) NSVfl (3) 
iKiminki et (12 00911: ( 41 RAITS (our data set 2); (5) iKiminld et 
(l200ii: fh'llSouzaet alJ (l2014h: tV l Miminki etal.l ll2012h : f8i lOMCl'’!: 
(91 ISazonov I lll961h : (10 1 iHalll (fT^4ll: (11) HIPPARCO^ (12) 


iLinder et alJ diod^ tl . 31 iMiczaikj 1 19531): (141 i Bohannan &Contil 


(119761): 1 ISllRauw et all fl999ll: (IhllWson & AbliTTgSlE 


“ http://skydot.lanl.gov/nsvs/nsvs.php 
https://sdc.cah.inta-csic.es/omc/index.jsp 
http://www.rssd.esa.int/SA/HIPPARCOS/apps/PlotCurve.html 


The Cygnus OB2 Radial Velocity Sur vey (RVS) of 1 28 se- 
lected O and early-B stars is underway d^minki et al.l l2007l 
I 2 OO 8 II 2 OO 9 I l2012t iKobulnickv et al.ll2014ll to assess the binary 
properties of massive stars within a common formation environ¬ 
ment at a similar age. There are now 48 massive binary systems 
known in Cygnus OB2. 

We have chosen Cygnus OB2 to study the mass loss in 
massive binaries because of the previously published data for 
massive stars in that association. Based on the radial-velocity 
(RV) data and our own photometric data that we obtained dur¬ 
ing 2011-2014, we have chosen seven massive binary systems 
with eclipsing light curves. In order to assess the change of pe¬ 
riod in those systems, we use a Bayesian parameter estimation 
method. We also use binary modelling software to acquire orbital 
and physical parameters. For several stars that were previously 
modelled using RV data only, we present models with additional 
photometric data. 


2. Data 

In this work we combine new photometric observations with 
photometry and radial-velocity data from the literature. 

References to the data that we used for each star are given in 
Table [1] We use the same naming conven tions as in the Cygnus 
OB2RVS papers: th e MT names fr om iMassev &^hompsonl 
(Il99lh. Schulte from ISchuit3 (Il958h and A36 and B17 from 
IComeron et al.l (l2002h . The data sets are described in more detail 
below. 

We monitored Cygnus OB2 with two telescopes during 
2011-2014. The first data set was obtained with the commer¬ 
cial iTelescope.net 0.25-m Takahashi Epsilon telescope (T4), lo¬ 
cated in Mayhill, New Mexico, USA. An SBIG ST-10 CCD cam¬ 
era with Johnson-Cousins BVI filters was used, providing a field 
of view (FOV) of 60.5 x 40.8 arcmin. Observations were car¬ 
ried out on 344 nights during 2011-2013. An observation con¬ 
sisted of three consecutive exposures in each filter. D ata reduc¬ 
tion w as carried out using the PHOTWORK software (iTuvikenel 
I 2 OI 2 I 1 . Magnitudes were extracted by means of aperture pho¬ 
tometry, except for MT720 that required PSF photometry due to 
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Fig. 1. The full frame of Cygnus OB2 region from our first data set 
(iTelescope) with the selected binaries shown in red. The dashed line 
indicates the FOV of our second data set (RAITS). 



a neighbouring star. To minimise observational errors, we used 
the nightly mean instrumen tal differential photometry values (for 
details see lFaur et al.ll201^ . 

The second data set was obtained on 103 nights during 
two seasons, 2013-2014, using the Tartu Observatory 0.31-m 
RAITS telescope, located in Toravere, Estonia. The telescope 
was equipped with an Apogee Alta U42 CCD camera and 
Johnson-Cousins BVRI filters. The FOV was 37 x 37 arcmin. 
Five consecutive exposures were taken with each filter during 
an observation. The CCD frames were reduced and magnitudes 
were extracted by aperture photometry as described above. 

Figure [T] shows the FOV of the two telescopes we used. Two 
of the seven program stars, A36 and B17, fell outside of the 
RAITS FOV. 

The Cygnus OB2 RVS is a spectroscopic survey of unbi- 
asedly chosen 128 stars of which 83 are B-type and 45 are O- 
type stars. They have used Keck-i-HIRES spectrograph (1999- 
2001), the Tick-(-Hamilton echelle spectrograph (1999-2000), 
the WIYN-i-Hydra spectrograph (2001-2008) and a longslit opti¬ 
cal spectrograph with a 2.3-m telescope at the Wyoming Infrared 
Observatory (2011-2013, 2014 May). The Cygnus OB2 RVS 
data that we used, ranges from 1999 to 2009. The longest data 
set of 43 measurements (1999-2007) is for MT059. Schulte 3 
and A36 have the least observations (~15) and a time base of 
one year while MT720 have 32 observations in two years and 
MT696 have 31 measurements with the time base of 30 days. 

The Northern Sky Variability Survey (NSVS) covers the en¬ 
tire northern sky in an unfiltered optical region, ranging from 8 to 
15.5 magnitudes fro m the Robotic Optical Transient Search Ex¬ 
periment (ROTSE-I) dWozniak et al.lbooil) . Except for MT696, 
there are 170 data points obtained within one year (April - De¬ 
cember, 1999) for all of our selected stars. The data set for 
MT720 was available but due to large photometric errors, the 
signal of its eclipses was unattainable and we had to exclude it 
from our analysis. The errors for NSVS data sets are typically 
larger than our own data sets because their data quality is limited 
by a severe blending due to 14-arcsec pixel size. 

The Hipparcos catalogue data (Ivan Feeuwenl200^ are avail¬ 
able only for the brightest star in our selection, Schulte 5, be¬ 
cause of the limiting magnitude of Hipparcos observations. The 
available data span over three years, from the end of 1989 to the 
beginning of 1993, with hundred broad-band measurements. 
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Fig. 2. The available data shown in time domain with reduced Julian 
dates at the bottom and calendar years at the top. Point sizes are indi¬ 
cating the number of measurements in a one-day bin. We have changed 
the scale of the plot in the middle (dashed line) to better show the data 
available from 1990 onwards. Green points indicate our photometric 
data sets, red points are the availa ble radial - veloci ty data, purple points 
show the times of minimum from ISazonovI l ll961h and black points are 
all other available photometric data listed in Table[T] 


The Optical Monitoring Camera (OMC) on board of the 
ESA International Gamma-Ray Astrophysics Laboratory (IN¬ 
TEGRAL) mission consists of 5°x 5° EOV photometric ca mera 
observing in the Johnson V-band (iMas-Hesse et alj[2003h . The 
only stars in our selection that have been observed by lOMC are 
Schulte 3 and Schulte 5. There have been more than 5000 obser¬ 
vations made between 2002 and 2009 (which is our longest time 
series from a single source) and their optical accuracy is roughly 
on the same level with NS VS photometry. 

ISouza et"^ (l2014l) have been monitoring MT696 since 2010 
in two 5-nm wide filters centred on the continuum and Ha (645 
and 656 nm respectively). Their photometric data span for over 
three years from 2010 to th e end of 2013 and co nsist of 105 
observations. These data and iKiminki et al.l (1201 2h RV data are 
the only published time series of MT696 to date. 

Eor Schulte 5 we have the longest time base with measure¬ 
ments from 12 different sources. The furthest published mea¬ 
surements for Schulte 5 are the times of minimum from 1898 al¬ 
though the accuracy of th ose measureme nts are under suspicion 
by the author s themselv es ( Sazonovll961h . During a 5-month pe¬ 
riod in 1952. lMiczaikal(11953 ) obtained 36 measurements with a 
photoelectric equipment o n the 32-cm refractor at the Heidelberg 
Observatory. iHalll (Il974l) observed Schulte 5 photoelectrically 
for 142 times with a 60-cm Seyfert reflector during the 1967 to 
191 0 period using a set of closely matching UBV passband fil¬ 
ters. fLinder et al.l (l2009h acquirecl over 300 exposures during a 
14-night campaign in 1998 with five filters of which three were 
centred on different emission-line wavelengths found in various 
subsets of Wolf-Rayet objects and other two were centred on 


line-free areas (around 5057 and 6051 A). Eor our period change 
analysis we chose only the measurements taken with the contin¬ 
uum narrowband filters as the emission line data also indicate 
variability not directly related to the eclipses of the binary. 

Radial-velo city data for Schu lt e 5 are available from 
three sources dWilson & A^ 1 195 it iBohannan & Conti! Il976t 
iRauw et al.lll999h . All the observations span three years (1946- 
1949, 1971-1974 and 1994-1997) and consist of 42, 7 and 25 
measurements accordingly. Schulte 5 has the longest time base 
for RV data as it is the only binary in our selection that has had 
more than one RV observation campaign carried out. 



Fig. 3. The V-filter phase diagram of Schulte 5, fold e d with a constant 
period. The iTelescope data (black points) and iHalll (11974) data (blue 
points) have been obtained 45 years apart. Solid lines are smoothed em¬ 
pirical light curves. The visible phase shift is due to the cumulative ef¬ 
fect of period change {P). 


All the available data points are shown in Fig. |2] to illus¬ 
trate the time base for each star. It can be seen from the figure 
that Schulte 5 is by far the most observed star in our selection 
whereas the other massive binaries have been scarcely studied. 
Our analysis depends mostly on three sources: our own two data 
sets from 2011 to 2014, NSVS and Cygnus OB2 RVS. While 
Schulte 5 has a time base of 116 years, all the other binaries 
have a shorter time base than 15 years. 

The time base is important in period change studies as a con¬ 
stant change in period introduces a gradual shift in eclipse tim¬ 
ings for each consecutive orbit. This causes a measurable cumu¬ 
lative phase shift over a long time base. The longer the time base, 
the more apparent the shift becomes, increasing the accuracy of 
the period change assessments. Figure |3] illustrates the cumula¬ 
tive effect of the period change for Schulte 5 with two different 
data sets that are observed 45 years apart. 


3. Method 


To estimate the variations in orbital period of our selected mas¬ 
sive binary systems, we used a Bayesian approach of simultane¬ 
ously fitting all available observational data: light curves, times 
of minimum and RV data. 

This method consists of two distinctive steps: acquiring the 
model curves and then assessing the period (R) and period 
change (R) values that simultaneously fit the data to those curves 
(Table EJ. We assume that the shape of the model curves re¬ 
main constant throughout the whole data serie^l For the ini¬ 
tial asse ssment of R, we us ed a phase dispersion minimisation 
method (IStellingwerflll97^ . This initial period was necessary 
to fold the photometric and RV data for input to the modelling 
software. For the modelling we used PHysics Of Eclipsing Bina- 
riEs (PHOEBE), (legacy version 0.31a), a bi nary modelling soft¬ 
ware based on the Wilson-Devinney code (IWilson & Devinn^ 
197lh. The modellin g was carried out using RV da t a from 


odellmg was carried out using KV da t a irom 
(120081 120091 12012h : IBohannan & Contil (Il976h : 


Kiminki et a . _ _ _ 

Rauw et al. ( 19991) and photometric data from our 2011-2013 


data set. We used this particular photometry data set as it had rel¬ 
atively smaller errors and included all of our s elected stars. The 
spectr al classes for the objects were taken from iKobulnickv et al.l 
(12014l) and they were used as an input for the initial temperature 


’ For more accurate modelling, apsidal motion should be taken into 
account. 
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(lGravll20()l!) and mass values (iMartins et al.ll2005h for the mod¬ 
els. The resulting model parameters for masses (Mi,M 2 ), mass 
ratio (q), radii {Ri,R 2 ), temperatures (Ti, T 2 ) and orbital param¬ 
eters are presented in Tabled 

The phased curves from PHOEBE software were used 
as the input to the Bayesian inference tool ca lled MultiNest 
dEeroz & Hobsonl[2008t lEeroz et alJ[2009l l2013h . The Bayesian 
parameter estimation method is a Monte Carlo algorithm that 
calculates the ephemerides of the star as 

H]D ^H]Do +PE+^PPE^, (1) 

where E is the cycle number. 

Using a random walk approach, we assessed the values of P 
and P that best set the data to the model curve. By comparing 
the folded light curve and folded RV curve that were calculated 
using the ephemerides to the model curves, we get the likelihood 
of the parameter fit. The global maximum of the likelihood pos¬ 
terior is the best ht of the parameters in conjunction with the 
model curves. We assumed the posterior to be Gaussian and thus 
the parameter ht errors are given as a standard deviation. To min¬ 
imise the errors, we hxed the value of HJDq using the minimum 
of the modelled PHOEBE light curve. 

Eor the period change analysis, we used one model light 
curve for all the photometric hlters as the htting of the 
ephemerides depends mainly on the time domain as opposed to 
the small differences in eclipse depths from different hlters. 

To test whether the PHOEBE model phased light curve 
causes any systematic errors in the period change determination, 
we constructed an empirical model light curves for our selected 
stars. All the results for P and P were obtained using both the 
PHOEBE model and the empirical model. The empirical phased 
light curve which was sigma clipped and smoothed, was derived 
from all the photometric measurements where the hrst estima¬ 
tion of the period change was already accounted for. As there 
were not enough RV data distributed over the full orbital phase, 
no empirical RV model curve could be accurately made. 

The errors for each photometric data set were obtained sep¬ 
arately as the standard deviation of the difference between the 
data points and the empirical light curve. The errors for RV data 
were taken fro m the respective p apers and the errors for the times 
of minimum in ISazonovI (1196 Ih were given as the timing error 
of the observations. These errors determined the weight of each 
data set in the hnal period change analysis. 

Conventionally, O - C diagrams are constructed using only 
precisely measured photometric times of the eclipses where the 
times are obtained by continuous monitoring of the selected star 
during an eclipse. Any RV data are usually completely ignored 
or only the times of conjunctions are used in the O - C anal¬ 
ysis. In our case though, RV data can be htted together with 
any photometric data. This method is especially useful if one 
has sparsely covered time series of the eclipses and does not 
know the exact times of the minima. This can happen either 
due to observational difficulties or when working with older data 
that have not been optimised in time domain for these kind of 
studies. In general, our method can be used for any data com¬ 
ing from all-sky-covering campaigns. They usually provide un¬ 
evenly spaced observation epochs that make it harder to detect 
the times of minimum but on the other hand have a long time- 
series for a lot of objects. Our simultaneous htting of data with 
the appropriate weights provides a solution that is not depen¬ 
dent on the measurement errors of any single minimum and also 
improves the statistical signihcance of the period assessment as 
more data points are used. Previous works for hnding period 


change using a ll of the available d ata have been carried out for 
RRd ty pe stars (iPaparo et al.||19 98h a nd for selected eclipsing bi¬ 
naries J Wilson & Van Hammdl201^ . 

4. Results 

We present the estimated orbital ephemerides in Table |2] 
and binary parameters in Table [3] We give one-sigma 
likelihood-estimation uncertainties for P and P, formal errors 
from PHOEBE for 7’2,f, e and calculated formal errorfl for 
M\,M 2 ,R\,R 2 and a. 

Our approach in this work does not take into account any 
cyclic patterns of period changes which may come from the light 
travel time effect through the presence of a third bo dy in the sys¬ 
tem o r from the stellar magnetic actiyity cycles (iLiao & OianI 
I 2 OIOI 1 . As data for most of our selected stars haye a relatiyely 
short time base compared to the typical periods of the third or¬ 
biting body or magnetic cycles, we haye opted to not include a 
cyclic pattern as a free parameter because this would artihcially 
introduce a new degeneracy to the fit. So any cyclic patterns in 
the period change are sources of possible errors in our results. 

Also, we haye not taken into account the apsidal motion in 
binaries with eccentric orbits (MT720 and MT059). The apsidal 
motion would yisibly alter the shape of light and RV curyes but 
we did not see this in our data sets (see Appendix lAl). 

Eigure 0] shows the parameter P and P estimation results 
for both empirical and PHOEBE model folded light curyes. 
The shape of the resulting parameter estimation likelihood cloud 
shows the amount of degeneracy of the two parameters as well 
as the uncertainties in the parameter estimations. Due to the de¬ 
generacy of P and P, the giyen uncertainties for those parameters 
are determined by marginalised distributions. If we had also fit¬ 
ted the HJDo yalue, the results would haye been degenerated be¬ 
tween all three of those parameters. Therefore, the yalue of HJDq 
was fixed to reduce the systematic errors of the two other htted 
parameters that we are interested in. The distance between the lo¬ 
cations of the parameter estimation likelihood clouds shows the 
systematic errors that arise from different model curye shapes 
and the sizes of those clouds show the statistical errors. So one 
has to pay attention to the model folded light curyes as they in¬ 
troduce some systematic errors compared to the true light curye. 

Appendix lAl describes the model hts for each of the dataset 
separately. The yisual comparison between the data and the mod¬ 
els can be used to assess the applicability of the final resulting 
parameters. 

By comparing the results to each other from the empiri¬ 
cal and PHOEBE folded light curye models for each star in 
Eig.m the results look consistent for most cases. The results for 
Schulte 5 seem to disagree with each other but, due to small 
uncertainties in the parameter estimation, are actually in good 
agreement. The resulting uncertainties for Schulte 5 are yery 
small because its time base is the longest out of our selected 
objects. The results for MT720 show a yery large systematic er¬ 
ror between the two models. This is due to the relatiyely short 
time base as well as the difference between the shape of those 
two model folded light curyes (see also Eig. lA.4l) . 

Eor four stars, we detect a non-zero period change yalues. 
The strongest eyidence for period change is found for Schulte 5, 
while MT696 also shows period change at three-sigma confi¬ 
dence leyel and A36 and MT059 at two-sigma confidence leyel. 
MT059 and MT696 exhibit a negatiye period change, while A36 
and Schulte 5 haye a positiye period change. Our calculations for 

^ http://phoebe-project. 0 rg/l.O/docs/phoebe_manual.pdf 
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Table 2. Orbital ephemerides from the Bayesian parameter estimation algorithm. 


Parameter 

A36 

B17 

MT059 

MT696 

MT720 

Schulte 3 

Schulte 5 

HJDo 

Temp (d) 

Pemp (X lO-**) 
Pmodel (d) 

Pmodel (X 10“*) 

56282.520 

4.676040(21) 

3.8 ± 1.3 
4.676040(19) 

3.9 ± 1.3 

56383.961 
4.021767(38) 
1.1 ± 1.6 
4.021719(39) 
-0.8 ± 1.7 

56535.672 
4.852281(39) 
-4.7 ± 1.7 
4.852265(26) 
-3.6 ± 1.2 

56181.719 
1.46918869(22) 
-0.37 ± 0.08 
1.46918857(24) 
-0.46 ± 0.08 

56585.602 
4.3619709(98) 
8.1 ±2.1 
4.3617773(96) 
-11.7 ± 3.2 

56264.633 
4.745964(15) 
-1.5 ± 1.0 
4.745998(29) 
0.0 ± 0.9 

56182.906 
6.597989(24) 
1.16 ±0.03 
6.5980120(30) 
1.30 ± 0.03 


Notes. Listed in this table are the epoch of periastron (HJDo), period (P) and the rate of period change (P). The parameters P and P are calculated 
using empirical (emp) and PHOEBE (model) phased light curves. HJDq is fixed in both cases. 
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Fig. 4. The parameter estimation results for period and period change. 
Periods are given as a difference from the Pe,^p value (see Table El- 
Coloured solid line represents the results for an empirical phase model 
and dashed line represents the results for calculated PHOEBE phase 
model. Inner contour is the 68% and outer contour is the 95% confi¬ 
dence level. 


B17, MT720 and Schulte 3 do not indicate any change in their 
periods. 

The stellar parameters (Table O are obtained with the 
PHOEBE software where RV data are fitted alongside photo¬ 
metric data. The stellar parameter values are not fine-tuned for 
each individual star as they are a by-product of acquiring model 
phase curves. As the parameter modelling is degenerate between 
many of the parameters, adding high-resolution spectra for pre¬ 
cise temperature and spectral-type measurements would con¬ 
strain the models further. 

There are no parameters for B17 because the only available 
RV data from IComeron et al.l (l2002h are published as a phased 
graph with no available ephemerides. Without RV data, no or¬ 


bital parameters or masses can be accurately modelled. Also, for 
MT059, the RV data are only available for the primary compo¬ 
nent and thus the model stellar parameters are only one of the 
possible parameter fits due to the degeneracy of the mass ratio. 
All the fits of our data to the acquired binary models are pre¬ 
sented in Eigs. m-M 


5. Discussion 

We use our estimated period change values to derive mass-loss 
rates for two distinctive cases. If we assume a conservative mass 
exchange between the two components with constant system or¬ 
bital angular momentum and system mass, the mass tra nsfer rate 
can be calculated as case I in ISingh & Chaub^(ll986h : 


1 M1M2 P 

3 Ml -M 2 P' 


( 2 ) 


We can also ass ume that the mass loss is solely through stel¬ 
lar wind (case IV in ISingh & Chaub^ll986l) : 


Ml =-i(Mi+M2)|. (3) 

Table |4] summarises for all of our selected stars the primary 
component mass-loss rates calculated for both conservative mass 
exchange and mass loss through stellar wind using the two ac¬ 
quired period change values of Pemp and Pmodei- Equations (|2]) 
and (O are for the idealised occurrences where no other effect 
interferes with the period-change value. The primary component 
mass-loss values in Table |4] represent isolated cases through ei¬ 
ther conservative mass exchange or stellar wind. They do not in¬ 
dicate the real values of the mass loss as in reality the mass loss 
is a combination of different factors. To constrain the mass loss 
values further, some additional method to find the strength of 
stellar wind (radio observations or analysis of the stellar-wind- 
induced spectral line profiles) has to be included. 

In the next sections we compare our results to the previous 
findings as well as present some of the problems that we came 
across on a star by star case. 


5.1. Schulte 5 


Schulte 5 (V729 Cyg, B D -h 40° 4220) is an eclip sing binary of 
two O-type supergiants dKobulnickv et al.ll2014l) . There is also 
evidence for a third non-the rmal radio source around the binary 
with a period of 6.7 years ([Kennedy et alJl20TO : ICazorla et al.l 
12014. 


Linder et al. (l2009t) assumed the orbital solution from 


iRauw et 4 ( 19991) to get the absolute masses of the components 
(Ml =31.9 + 3.2 Mq and M 2 = 9.6 + 1.1 Mq). Our modelling 
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Table 3. Binary parameters from the PHOEBE modelling software. 


Parameter 

A36 

B17 

MT059 

MT696 

MT720 

Schulte 3 

Schulte 5 

M, (Me) 

22.41 ± 2.96 

- 

21.23 ± 6.24 

14.00 ± 0.81 

18.52 ± 1.51 

32.73 ± 7.04 

38.18 ±5.21 

M 2 (Mo) 

15.32 ± 2.02 

- 

6.22 ± 1.83 

12.67 ± 0.73 

13.19 ± 1.07 

17.21 ± 3.70 

10.24 ± 1.40 

q 

0.68 ± 0.04 

- 

0.29 ± 0.01 

0.91 ± 0.03 

0.71 ± 0.03 

0.53 ± 0.05 

0.27 ± 0.03 

Rl (Re) 

16.08 ± 0.27 

- 

12.51 ±0.10 

5.92 ± 0.09 

9.93 ± 0.59 

18.40 ± 0.41 

27.16 ±0.34 

R 2 (Re) 

13.48 ± 0.30 

- 

5.89 ±0.21 

5.63 ± 0.09 

8.47 ± 0.38 

13.50 ± 0.34 

15.15 ±0.20 

Ti (K) (fixed) 

27000 

- 

38000 

32000 

29800 

41200 

36100 

72 (K) 

22924 ± 38 

- 

18966 ±98 

31070 ± 25 

19202 ± 80 

29702 ± 37 

23681 ± 112 

if) 

71.76 ±0.11 

- 

68.18 ±0.13 

81.87 ±0.06 

71.09 ±0.28 

57.47 ± 0.09 

64.65 ± 0.05 

a (Re) 

39.43 ± 1.45 

- 

36.34 ± 3.46 

16.23 ± 0.24 

35.52 ± 0.78 

43.72 ± 2.63 

53.90 ± 2.08 

e 

0.01 ± 0.001 

- 

0.14 ±0.002 

0 (fixed) 

0.34 ± 0.005 

0.015 ± 0.001 

0(fixed) 

w (°) 

350 

- 

189 

0 (fixed) 

298 

0 

0(fixed) 

Spectral class 

B0Ib±B0III 

07: ±09: 

08V±B 

O9.5V±B0V 

B0-B1V±B1-B2V 

06IV: ±09111 

07Ianfp±0fpe/WN9 


Notes. Parameters in this table are the primary and the secondary component masses (Mi, M 2 ), the mass ratio M 1 /M 2 (q), the primary and the 
secondary component radii (Ri, R 2 ), the primary and the secondary component temperatures (Ti, T 2 ), the inclination ang le (;), the semi-ma j or axi s 
(a), the eccentricity of the system (e), the argument of periastron (tu) and the binary s pectral clas ses that are taken from iKobulnickv et al.l l l2014h . 
The temperature for the primary component was fixed according to the spectral class llGravll20()Ih . 


Table 4. The mass transfer and mass loss through stellar wind rates for 
our selected stars. No mass loss through stellar wind is given for the 
negative P value as this would indicate a positive mass loss from the 
primary. The errors are calculated using the P, P, Mi and M 2 uncer¬ 
tainties. The errors for MT696 are large due to the similar companion 
masses. 


Star 

used P 
value 

Ml mass transfer 
(IQ-^Meyr-') 

Ml stellar wind 
(10“®Moyr-') 

A36 

p 

^ emp 

48 ± 19 

-56 ± 19 


P model 

50 ±21 

-58 ± 19 

MT059 

p 

^ emp 

-10.5 ±4 

- 


P model 

-8.0 ± 3 

- 

MT696 

p 

^ emp 

-41 ± 160 

- 


P model 

-50 ± 254 

- 

MT720 

p 

^ emp 

104 ± 30 

-10.7 ±4.2 


P model 

-149 ±45 

- 

Schulte 3 

p 

^ emp 

-13.9 ± 10.2 

- 


P model 

0.001 ± 8.8 

- 

Schulte 5 

p 

^ emp 

3.0 ± 0.2 

-15.6 ±0.7 


P model 

3.4 ± 0.2 

-17.4 ±0.8 


Notes. The two period change values (Pemp and Pmodei) for each star are 
given in Table |2] B17 is omitted from the Table as there are no mass 
estimates. 


reveals similar masses for the components, 38.2 + 5.2 Mq and 
10.2 + 1.4 M 0 accordingly. 

iMartins et al.l (l2005h hav e calculated a grid of m odels for O 
type stars using CMFGEN (iHillier & Millel^f^998l) atmosphere 
code and they include theoretical spectroscopic masses for them. 
For the primary component of spectral type 071, the spectro¬ 
scopic mass would be Mjpec = 40.91 Mq. They note that one 
should take care in dealing with their spectroscopic modelled 
masses as the uncertainty may be as high as 35 to 50%. Never¬ 
theless, our resulting mass for the pri mary is more close t o the 
theoretical value than to the result from iLinder et alJ (l2009l) . 

We chose Schulte 5 as a testbed for our period change study 
as the eclipse timing method for period change as well as the 
mass-loss rates have been studied on multipl e occasions. Previ¬ 
ous assessments for the mass-loss rates are bv iPersi et al.l (Il990b 
derived from infrared ob servations M - 2.5 x 10“^ Moyr^', by 
IConti & HowarthI d 19991) combined with multi frequency radio 
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Phase 


Fig. 5. The binary modelling for Schulte 5 with the PHOEBE software. 
Upper panel shows our photometric iTelescope data (data set 1) with 
the model curve (s olid line) and lower panel sho ws combined radial - 
velocity data from [Bohannan & Conti ( 119761) and iRauw et 
with the model curves for both primary and secondary components. 

observ ations M - (3.7 + 1.3) x 10“^ Moyr^' and bv ETnder et al.l 
(l2009h using the period change rate of observed times of mini¬ 
mum M= (2.1 + 0.6) X 10“^ Mgyr"*. 

Fo llowing the calculation from Equation (|2|i, iFinder et al.l 
(l2009h found the mass-loss rate Mi = (4.7 + 1.5) x 10“®Moyr“' 
which they argued against as this value was lower than ex¬ 
pected from the previous assessments derived from the infrared 
and radio data. Using the masses obtained from our PHOEBE 
model and the P and P from our Bayesian parameter esti¬ 
mation method, we also get similar mass-loss rates of M\ - 
(3.0 + 0.2) X 10^^ Moyr^‘ and Mi = (3.4 + 0.2) x 10^® Moyr^' 
dependit ig on the applied p eriod change value. This difference 
between [Linder et al.l (l2009l) Mi value and our values is mostly 
due to the differenc e in the period change values. The value from 
iLinder et'aD (l2009l) . P = (1.9 + 0.5) x 10“*, is higher than our 
obtained value of Pemp = (1-16 + 0.03) x 10“* and thus their 
mass-loss rate is also a little bit higher. Our period change value 
should be more accurate as we have over 40 years longer time 
base. 
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Fig. 6. The binary modelling for A36 with the PHOEBE software. Des¬ 
ignations of lines, labels and photometric data set are the sam e as in 
Fig .|5] Radial-velocity data are taken from iKiminki et all ll2009h . 


lYasarsov & Yakut! (I20l4) found a period change of 
PIP ~ 9.4 X 10^ from eclipse timings by acquiring new ob¬ 
servations on 65 observing nights from 2010 to 2011 in Bessel 
L® VK/bands. With their given period of P = 6.597981d, we can 
calculate P - 7.01913 x 10“^, which is much higher from the 
value of P = (1.16 + 0.03) x 10“*^ that we obtained. It may be 
possible that their given P/P value is in the units of years instead 
of days as one might assume from their period value. In this case 
the value P = = 1.92 305 x 10“® agrees with the value of P from 
iLinder et alJ (l2009t) . They also model the masses of the primary 
(31.6 + 2.9 Mo) and secondary (8.8 + 0.3 Mq) component. If 
we use their data for calculating mass loss through mass transfer, 
we get Ml = 4.32 5 x 10~^ Moyr~^ w hich is very similar to the 
mass-loss rate that iLinder et al.l (l2009t) obtained but differs from 
ours due to the value of period change and obtained component 
masses. 

We find that the mass-loss rate through stellar wind (Eq. Q) 
for Schulte 5 is Mi = (-1.56 + 0.07) x 10“^ Moyr"' or 
Ml = (-1.74 + 0.08) X 10“^ Moyr“' depending on the P 
value, which are comparable to the result of iLinder et alJ (l2009ll 
(Ml =(-2.1 + 0.6) X 10“^ Moyr^'). The previous mass-loss 
rate assessments of the stellar wind from radio and infrared ob¬ 
servations have confirmed the value to be close to our finding. 
This indicates the period change of the system to be caused by 
the mass loss from stellar wind of the primary component instead 
of the mass transfer from the secondary component. 


5.2. A36 


A36 was first named in a paper bv IComeron et al.l (l2002h where 
they listed it as a possible bi nary. The ec l ipsing nature of A36 
was found from NSVS data (lOterol120081: iHoffman et all 20091) 
and it is classified as an Algol/yS Lyrae candidate. Kiminki et al.l 
determine the spectral type as a BOIb-i-BOIII system on 


the basis that the light curve and spectra show similar component 
temperatures and luminosities with differences in specific lines. 

Our calculated mass rat io for the binyy syst em is 0.66 and 
it is within error limits of iKiminki et al.l (120091) estimation of 
0.70 + 0.06. They give a period of P = 4.6752 + 0.0007 d that 
is slightly different than our result of P = 4.67604 + 0.00002 d. 


This is mainly due to the lack of observations in the RVS 
for A36 to precisely determine the period. We find the pe¬ 
riod change from the combination of two models to be P = 
(3.85 + 1.3) X 10“*, thus the mass loss by direct mass transfer 
would be Ml = (4.9 + 2.0) x 10“^ Moyr“' and by stellar wind 
Ml =(-5.7 + 1.9) X 10“5Moyr“'. 

The value of Mi by stellar wind is too high for a BOIb star 
as the maximum empirical est imates for mass-loss rates are typi - 
cally lower than 10“^ Moyr * (IVink et alJ2000t|Puls et alj200^ . 
Therefore, the change in period must be caused by the mass 
transfer from the secondary to the primary. This also indicates 
that the secondary component has to be more ev olved than the 
primar y. Although the spectral classification by IKiminki et al.l 
(l2009l) contradicts this conclusion, it may still be viable if the 
luminosity class of the primary component could be lowered. 
Our PHOEBE modelling results for the stellar radii and masses 
also support the argument that the primary component is more 
close to the luminosity class III. 


5.3. B17 

B17 was first discussed in a paper bv IComeron et al.l (l2002h 
where they discovered that it displays emission in the Brackett y 
line (the B notation in front of the star number) which is an indi¬ 
cator of high mass-loss rate, typical to the evolved stellar stages. 
It is a massive evolved binary system wit h a spectral cla s sifica- 
tion estimation of 07Ia-i-09I supergiants. [Hoffman et al.l (l2009l) 
classified it as an Algol/j0 Lyrae type using NSVS data. 

Eor B17 we could not use any spectroscopic data. 
IStroud et al.l (l2010l) have done an extensive photometric and 
spectroscopic campaign on B17, but unfortunately their data are 
not publicly available and they are only in a phased form on a 
graph. Their determined period of 4.02174 + 0.00003 d coin¬ 
cides with ours and as the time base between the NSVS pho¬ 
tometry data and their observed spectra was only four years, 
they do not discuss any possibility of a period change. We 
did not find any indication of a period change within the bi¬ 
nary system with values of Eemp = (1.1 ± 1.6) x 10“* and 
Emodei = (-0.8 + 1.7) X 10“* and therefore we find no mass 
loss in the system. Also due to no RV data available, no stellar 
parameters were modelled. 


5.4. MT059 

MT059 (Schulte 1) is OS V type star wi th a B, possibly later 
than BIV, type companion. IKiminki et al.l(l2008l) have only mea¬ 
sured radial velocities for the primary component and have not 
detected any spectral features of the companion thus indicating 
that the luminosity ratio has to be less than L 2 /T 1 < 0.2. Eor 
their primary component mass assess ment, they use t he the oreti- 
cal mass for a spectral type 08V from iMartins et al.l (l2005h with 
a value of Mjpec = 21.95 Mq. They come to the conclusion that 
due to the absence of spectral signatures of the secondary com¬ 
ponent, the interpolation from its spectral type yields the mass 
of the secondary companion to be 5.1 Mq < M 2 < 13.8 Mq 
and the mass ratio 0.24 + 0.02 < q < 0.64 ± 0.03. 
Our resulting mass ratio falls on a lower side of their prediction 
with a value of 0.29 and the masses we obtained are in agree¬ 
ment with the previously mentioned predictions with values of 
Ml = 21.23 + 6.24 Mq and M2 = 6.22 + 1.83 Mq. As there 
are no RV data for the secondary companion, we cannot assess 
the mass ratio independently and thus the mass ratio is degener¬ 
ate with the other stellar parameters. Also the semi-major axis of 
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Fig. 7. The binary modelling for MT059 with the PHOEBE soft- 
ware.Designations of lines, labels and photomet ric data set are t he sam e 
as in Eig.|5] Radial-velocity data are taken from iKiminki et aO ll2008h . 



Fig. 8. The binary modelling for MT696 with the PHOEBE software. 
Designations of lines, labels and photomet ric data set are t he sam e as in 
Fig .[5] Radial-velocity data are taken from iKiminki et ai](l2012h . 


the system yields large formal errors which lead to large errors in 
the mass assessments. The modelled light curve (upper panel in 
Fig.S indicates the presence of a reflection effect (phase ~ 0.2) 
that is only loosely supported by the data. A more thorough in¬ 
vestigation is needed to determine whether this effect is from the 
observational uncertainties or from the underlying physics of the 
system. 

MT059 is indicated bv iKiminki et al.l (l2008l) as a SB 1-type 
binary with no mentioning of its eclipsing nature. Our observa¬ 
tions, as well as NSVS photometric data, show clear eclipses 
from the photometric folded light curve. The binary has an ec¬ 
centricity of 0.14 which is also noticeable from the uneven tim¬ 
ings in phase space between the primary and secondary eclipse. 

We And that out of all of our selected binary systems, MT059 
has the highest period change value (Pemp = (-4.7 + 1.7) x 
or Pmodei = (-3.6 + 1.7) X 10“^). This is due to the RV ob¬ 
servations that are also visibly shifted in our PHOEBE mod¬ 
elling (Eig. |7]l. The mass loss through conservative mass ex¬ 
change with the results from our empirical model is Mi = 
(-10.5 + 4) X 10~® Moyr^'. As the period change is neg¬ 
ative, the mass-loss from stellar wind cannot be calculated us¬ 
ing Eq. Q because it would result in a positive Mi value and 
therefore the period change must be from mass transfer from the 
primary to the secondary component. 

5.5. MT696 

MT696 is a binary system of an 09.5V spectral type prim a ry and 
an early B type companion. I Rios & DeGioia-EastwoodI (l2004h 
were the first to report the eclipsing binary status of MT696, 
indicating t hat it might be a W U Ma type binary with a 1.46- 
day period. iKiminki et al.l (l2008l) argue that the early spectral 
types of the components, the absence of emission features, and 
a period over 1 day suggest MT696 to be /? Lyr type instead. 
They also improve the period estimation to 1.4692 + 0.0005 
days and estimate the inclination to be near 80° and the mass 
ratio q to be 0.85 + 0.03, which are close to our estimations. 
Their minimum mass estimations are 15.1 + 0.7 Mq for the 
primary and 12.8 + 0.5 Mq for the secondary component, 
while our modelling gives the total masses of 14.00 + 0.81 Mq 


and 12.67 + 0.73Mo, respectively. Erom iMartins et al.l (l2005l) 
the mass for spectral type 09.5V is Mjpec = 16.46 Mq. Our 
mass estimations are within errors of the minimum masses from 
iKiminki et al.l (l2008h but lower than the theoretical primary 
mass, so our resulting masses might be underestimated. An¬ 
other possibility is that the spectral types and thus the masses 
of the binary component s might not be accurately determined by 
iKobulnickv et al.l (l2014l) . 

ISouza et^ (1201 4h And the period to be 1.46919 + 0.00006 d 
and they did not And any period change in their two year 
long data set. We, on the other hand, And the period change 
at three-sigma confldence level with the value of Remp = 
(-0.37 ± 0.08) X 10-^ andP^odei = (-0.46 + 0.08) x lO^^ This 
is the smallest period change value amongst our studied binaries. 
Due to the small value for period change and the similar masses 
for the binary components, the errors for conservative mass-loss 
rates are very high, as the total error depends on the diflerence 
between the companion masses. The acquired mass-loss rate for 
direct mass transfer is then M] =(-0.4 + 1.6) x 10“"^ Moyr^' 
and it is dominated by the measurement errors of the stellar 
masses. 

5.6. MT720 

IKiminki et al.l (1201 2|) have observed MT720 and measured its 
radial velocities on 32 occasions, noting that no known photo¬ 
metric variability have been studied. They indicate the period 
of the system to be 4.3622 + 0.0003 days (compared to our 
= 4.36197 + 0.00001 and = 4.36177 ± 0.00001 
days) and an eccentricity of 0.35 ± 0.02, which is evident 
from the RV curve. They also And minimum masses from the 
RV curve and suggest the spectral type of the binary system 
to be BO-BIV + B1-B2V. Their calculated minimum masses 
of Ml > 15.5 Mq and M 2 > 11.1 M© are in accordance 

with our modelling results for the component total masses of 
Ml ^ 18.52 ± 1.51 M e and M 2 ^ 13.19 + 1.07 Mq. 

ISalas et ^ (1201 4l) have done a 1.5-year-long photometric 
survey in Cygnus OB2 and have included MT720 as a conflrmed 
eclipsing binary with a period of 4.3619 ± 0.0001 days. Our 
photometric data also show clear evidence of eclipsing nature 
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Fig. 9. The binary modelling for MT720 with the PHOEBE software. 
Designations of lines, labels and photomet ric data set are t he sam e as in 
Fig .|5] Radial-velocity data are taken from iKiminki et all ll2008h . 


of the binary. At first we note that in phase space the primary 
eclipse coincides with the RV curve when the more massive star 
is eclipsing the secondary component, indicating that it must 
be less luminous. After modelling the light curve, it turns out 
that because of the high eccentricity of the system, the primary 
eclipse depth lessens into what we see as a secondary eclipse. 
This effect is due to the combination of the high eccentricity and 
a low inclination angle resulting in a partial eclipsing of the pri¬ 
mary component by the secondary star. The model in Fig. |9] is 
acquired by introducing a phase shift of 0.45 so the secondary 
eclipse would coincide with a phase zero point. 

Modelled stellar parameters for MT720 are more prone to 
errors compared to the other binary systems due to the bad fit 
between the data and the model. As can be seen in Fig.|2 the sec¬ 
ondary eclipse in the observations is narrower compared to the 
model curve. Our period change modelling for this system have 
a strong systematic effect due to the differences between the two 
phased light curve models. The above mentioned modelling dis¬ 
cord is due to the fact that MT720 have the largest photometric 
uncertainties as it is the faintest object in our selection. Because 
of those systematic errors between the models, we cannot find 
any reliable evidence for the period change in MT720. 

5.7. Schulte 3 

Schulte 3 is a massive binary s ystem consisting of 0 6IV: and 
09III stellar type components. iKiminki et al.l (l2008h found a 
mass ratio of q = 0.44 + 0.08 which agrees within errors with 
our result of 0.53 + 0.05 but differs from the th eoretical mass 
ratio o f q = 0.64 for 06IV; + 09III binary system (iMartins et alJ 
l2005h . We find our modelled masses (Mi = 32.73 ± 7.04 Mq 
and M 2 = 17.21 + 3.70 Mq) to agree with their minimum mass 
predictions of Ml > 17.2 Mq and M 2 > 7.6 M©. Our acquired 
primary component mass is also close to the theoretical spec¬ 
troscopic mass of 06V (Mspec = 31.73 Mq), but the theoretical 
mass value for the secondary component (Mspec = 23.07 Mq) is 
higher than our result. 

IKiminki et alJ (l2008b argued that the difference between the 
mass ratios is due to the mass transfer between the compo¬ 
nents. The period we find coincides with their result {P - 
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Fig. 10. The binary modelling for Schulte 3 with the PHOEBE software. 
Designations of lines, labels and photomet ric data set are t he sam e as in 
Fig .[ 5 ] Radial-velocity data are taken from iKiminki et all (l2008h . 

4.7459 + 0.0003 days) within the errors, but we do not detect any 
notable period change in the system. Our period change value for 
empirical model is Pemp = (-1.5 + 1.0) x 10“^, which does not 
imply a non-zero value at two-sigma confidence level, and for 
PHOEBE model the value is Pmodei = (0.0 + 0.9) x 10'*. These 
values show clear evidence for no period change within the sys¬ 
tem and indicate no mass transfer between the components. 

6. Conclusions 

We studied the period change of seven massive binaries in the 
Cygnus OB2 association. Eor this we used archival photome¬ 
try and RV data together with our own photometric observations 
from 2011-2014. Our observations were essential in the period- 
change analysis for A36, B17, MT696 and MT720 as they in¬ 
creased the time base up to 13 years. We note for the first time 
the eclipsing nature of MT059. 

Eour stars show evidence of a non-zero period change at 
two-sigma confidence level. For three stars we report the pe¬ 
riod change for the first time. Our calculations confirm the pre¬ 
vious assessments for P of a star Schulte 5. We also show that 
the change in period in Schulte 5 originates from the stellar 
wind from the primary component. Two of the systems, B17 and 
Schulte 3, show period change values close to zero indicating 
that there is no mass transfer between the binary components or 
that the period change values are too small for us to estimate 
them reliably with the available data. Our method had system¬ 
atic errors in fitting the period and period change parameters for 
MT720 and thus we did not derive any reliable values for that 
system. Two of the systems, MT059 and MT696, have negative 
P values that show mass transfer from the primary to the sec¬ 
ondary component. We argue that the change in period for A36 
might be from mass transfer from the secondary component to 
the primary, indicating a more evolved secondary stellar compo¬ 
nent. 

We show that applying the Bayesian parameter estimation 
method, period change can be calculated using all the avail¬ 
able data simultaneously for any binary star. Utilising all data 
points in contrast to using only the times of minimum, improves 
the statistical significance of the calculations and opens up new 
data sets to include in the analysis. This is especially important 


Article number, page 9 ofll2l 






























A&A proofs: manuscript no. Laur2015 


with the trend to digitise archives of photographic plates as it 
adds new historical observations that were not optimised for the 
eclipses of binaries, which means that the times of minimum 
may be difficult to obtain. 
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Appendix A: Individual data sets and model curves 

In this section we present the resulting phase diagrams for pho¬ 
tometric and RV data sets for our studied stars. The diagrams 
are constructed by folding the data as Equation ([T]i with the re¬ 
sulting orbital ephemerides (Table |2]l. Figures lA.l I - IA.81 show 
individual data sets together with the phased model curves. 
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Fig. A.l. Phase diagrams for A36 using the Pemp and Pemp from Ta¬ 
ble [2] Blue dashed lines in the magnitude plots denote a model light 
curve from PHOEBE software and the red solid lines denote an empir¬ 
ical light curve. The solid line in the RV plot represents the model for 
the primary component and the red dashed line represents the model for 
the secondary component. 
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Fig. A.2. The phase diagrams for B17. Designations are the same as in 

Fig.lXn 
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Fig. A.3. The phase diagrams for MT696. Designations are the same as 
in Fig. IA.il 
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Fig. A.4. The phase diagrams for MT720. Designations are the same as 
in Fig. IA.il 
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Fig. A.5. The phase diagrams for Schulte 5. Designations are the same 
as in Fig. IA.il 
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Fig. A.6. The 0-C diagram for Schulte 5. Points denote the time differ¬ 
ence between the observed times of minimum and the calculated times 
of minimum using the Pemp and Pei^p from Table |2] Dashed line shows 
the 0-C value of zero. 
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Fig. A.7. The phase diagrams for Schulte 3. Designations are the same 
as in Fie. IA.il 
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Fig. A.8. The phase diagrams for MT059. Designations are the same as 
in Fie. IA.il 
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